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Abstract 

Considering the nuclear motion, the authors give out the nonrelativistic ground energy of a 
helium atom by using a simple but effective variational wave function with a flexible parameter k. 
Based on this result, the relativistic and radiative corrections to the nonrelativistic Hamiltonian 
are discussed. The high precision value of the helium ground energy is evaluated to be —2.90338 
a.u., and the relative error is 0.00034%. 
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As is well known, studying the typical Coulomb three-body bound-state problem such 
as a helium atom is a very basic and important problem in the field of atomic physics 
ever since. Much work has been carried out by some approximate methods. One of the 
earliest variational calculations has been performed by Hylleraas [1] in 1929. Since this 
time, many physicists have contributed incremental steps in this endeavor, often trying 
to use mathematical insight for advantage j^-Q]. The strongest line of theory has been 
focused on the analytic properties of the wave function, especially following the work of Fock 
in 1954 2], which showed a weak logarithmic singularity at the three-particle coalescence. 
The recent works of Korobov and Charles Schwartz ect get very high precision result for the 
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But most of these authors mainly lay their works on the precision and convergence of 
the variational wave function except for the correction to Hamiltonian of a helium atom. 
Even the motion of nucleus is not considered, so the variational results are lower than the 
experimental value . In this paper, we first construct a simple but effective variational 

wave function with a flexible parameter k. Then, considering the motion of nucleus and the 
relativistic and radiative effects, we give out all the corrections to the nonrelativistic ground 
energy. Finally, the ground energy of a helium atom is worked out with the variational 
method and perturbation method. 

It is well known that the Schrodinger ground energy E and the corresponding wave 
function U are found as a solution of the variational problem 

f U*HUdr 

E = mm J rTT — . (1) 
u JU*Udr 

For a helium atom, we use a simple but effective variational wave function, which contains 
only even powers of t, as follows 

U( ri , r 2 ) = <p(ks, kt, ku) = e~ ks C lmn (ks) l (kt) 2m (kuT, (2) 

Iran 



s = n + r 2 , t = — n + r 2 , u— ri2, (3) 

where r x and r 2 denote the positions of the electrons with respect to the nucleus; s, t and u 
are called Hylleraas coordinates [lj; k is a flexible scaling parameter. The function may be 
expected to converge rapidly for medium-sized of s, t and u. The Hamiltonian is taken in 
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following nonrelativistic approximation in many papers 
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2 2 n r 2 ri 2 
where p x 2 are the momenta of the electrons and Z = 2 is the nucleus charge (in units of 
the proton charge). And the variational results in References Q]- ^ are lower than the 
experimental value, which shows that the Hamiltonian form in Eq. (J3J) is not accurate 
enough. So the movement of the nucleus should not be ignored, and the nonrelativistic 
Hamiltonian can be represented as follows 

P? p2 P 2 Z Z 1 

// » = Y + Y + 2M-n-^ + ^ (5) 

Here P = — p 1 — p 2 is the momentum of the nucleus and the finite nucleus-to-electron mass 



ratio M = m a /m e = 7294.299508(16). From the later computational results, it can be seen 
clearly that the variational energy is higher than the experimental value. 

To get more accurate result, we take into account the relativistic and radiative effects. 
The following operator describes the Breit a 2 correction to the nonrelativistic Hamiltonian 

Q 

H (2) = H 1 + H 2 + H 3 + H A + H s , (6) 
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H 1 = --(pt + pt), 
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Ho = — 



Pi • P 2 , (r 12 • PiXrig • p 2 ) 
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Ha = — (Pi Li + p 2 -L 2 ), 



L 2 x p 2 + 2 



r 2 i x p 1 
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(8) 
(9) 
(10) 



Hk 



a 



~(Si • S 2 )5(r 12 ) + i- 
3 r J 
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Si • s 2 



3(S 1 -r 12 )(S 2 -r 12 ) 
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Here Li = — ViV" (V = Z/r± + Z/r 2 — 1/7*12) is the Coulomb field due to nucleus plus the 
second electron. The physical significance of the various terms in Eq. (jSJ is as follows: 

Hi is the relativistic correction due to the "variation of mass with velocity" (which does 
not depend on electron spin). H 2 corresponds to the classical relativistic correction to the 



interaction between the electrons. This correction is due to the retardation of the electro- 
magnetic field produced by an electron. H3 is the interaction between the spin magnetic 
moment and the orbital magnetic moment of the electrons (spin-orbit coupling). if 4 is a 
term characteristic of the Dirac theory, which is also present in the Hamiltonian for a single 
electron in an electric field. H 5 represents the interaction between the spin magnetic dipole 
moments of the two electrons. 

For the ground state of helium-like atoms, I, s and j are all zero and there is no fine 
structure splitting. Nevertheless, the operators Hi to H 5 in Eq. © contribute relativistic 
corrections to the nonrelativistic energy eigenvalue of relative order {Za) 2 and Za 2 . These 
corrections are only the leading terms in an expansion in powers of a and Za. And we 
explicitly take into consideration that the spin of the nucleus and the total spin of electrons 
are both equal to zero. In particular, we replace the product of the electron spin operators 
Si • S 2 by its eigenvalue in the single state, —3/4. The expectation values of the operators 
H\ to H§ for the ground state of a helium atom with nuclear charge Z are 

a 2 

E 1 = --(<pt> + <pi>), (12) 



* 2 = -^ 



' r Pi P2 ^ „ ( r i2 ■ Pi)(ri2 ■ P 2 
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n 2 r» 2 



(13) 



E 3 = 0, (14) 



E, = ,a 2 <Z 6 -^M-5(r 12 )>, (15) 

E 5 = 2na 2 < 8{r 12 ) > . (16) 

The angle brackets in Eqs. (J 12)) to (|16|) and below denote the average value over the non- 
relativistic ground state variational wave function. It can be shown that the expectation 
value E2 of the operator H 2 vanishes if any wave function of product form U = u{ri)u{r2) is 
used, i.e. both for the hydrogen-like and for the Hartree wave functions. If a more accurate 
function, which includes the effects of polarization, is used, then a finite (but numerically 
small) value is obtained for E 2 . To simplify the presentation, we take E 2 = 0. So, the total 
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Breit a 2 correction to the nonrelativistic value is 

SE m = „ 2 < _p!±p| + rZ jW+jW + ri(ra) > . (17) 

8 2 

Order a 3 correction to the nonrelativistic energy can be represented as follows 19] 

SE m = a *I^f_2hi a -l3 + ^\ <6(tx)+ 6(t 2 ) > 

14 , i 1QA\ r/ , : 7 lnri 2 + 7 
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+ ( T In a + — 1 5(r 12 ) +-< — ^in • p > } . (18) 



Here 

n = ^, p = -iV = -i^, (19) 
ri 2 ot 12 

7 = 0.5772 is the Euler constant and /3 = 4.3700392 is the helium Bethe logarithm 20]. 

Now, we give out the calculations and results of the helium ground energy. 

Firstly, we shall discuss the ground state expectation values of the Dirac delta-function, 
the square of the kinetic energy operator and the nonrelativistic energy, which evaluated by 
using various wave functions. 

The expectations of the three-dimensional Dirac delta-function 5(ri) and <5(r 12 ) are 

00 

<S(t 1 )>= J dr 2 U 2 (0,v 2 ) = 4vr J r 2 (p 2 (r,r,r)dr, (20) 





< £( ri2 ) >= J dnU^TuT!) = An J rV(2r,0,0)dr. (21) 



It is very convenient to evaluate them in spherical coordinates. 

For the square of the kinetic energy operator for electron 1, or p| = Vf = Af, we first 
note the following equation 

Vi • [*7Vi(AiC/) - (A^ViU] = U(AjU) - (A^) 2 . (22) 

For any analytic function U that falls off exponentially at large distances, the integral of 
the left side of Eq. ()22|) over the whole ri-space must vanish (from Gauss' theorem). Then 
there are two alternative forms for the expectation value of pf 



< 



>= J dnd^UAjU = J dTxdT^AxU) 2 . (23) 
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Great care must be taken if the first form in Eq. ({23}) is used. For the exact wave function, 
AiU behaves like Zjr\ or like I/V12 if n or ri2approaches to zero, just like the potential 
energy in the total Hamiltonian. A\U then has a delta-function type of singularity at T\ = 
and r 12 = and a wrong answer would be obtained if the first integral in Eq. (J23|) were 
evaluated by a limiting process which excludes an infinitesimal region around the origin and 
around rxi = 0. The second form of Eq. (|23|) is free from these difficulties and is, in any case, 
easier to evaluate in practice. Considering the exchange symmetry of the wave function, we 
have < 5(rx) >=< 5(r 2 ) > and < pf >=< \>\ >. 

The expression of the nonrelativistic ground energy has the following form 

E =< H„ >= **» + £*, (24) 

where K is kinetic energy, P is potential energy and W is a normalization factor. The 
function Eo(k) is a quadratic parabola that take the minimum at k — —P/2K. When the 
number of basis functions N is large enough, the value of k will fluctuate between 2.0451 
and 2.0452. For example, we get k = 2.0451486913735 when N = 50. 

Secondly, we will present the computational values for all the contributions to the helium 
ground energy in atomic units (m p =l, h=l and e=l). The particular numerical results and 
their comparisons with the experimental value are collected in Table 1 and Table 2. E^ and 
Eq are calculated by variational method, 5E c h r , 5E^ and 5E^ by perturbation method. 



Table 1. Nonrelativistic energies E^ &ndE for the helium ground state and their 
comparisons with the experimental value. E^ corresponds to the expectation value of H^, 
Eq to the ex pec tation value of Hq and E exp to the experimental value, which equals to 



-2.90338629 



2l| (in atomic units). Ai^ = E Q 



E exp and AE 



En — E, 



exp- 



N 




AEm 


E 


AE 


20 


-2.90370938 


-0.00032309 


-2.90328962 


0.00009667 


30 


-2.90371945 


-0.00033316 


-2.90329969 


0.00008660 


40 


-2.90372103 


-0.00033474 


-2.90330128 


0.00008501 


50 


-2.90372124 


-0.00033495 


-2.90330389 


0.00008240 



Table 2. The relativistic and radiative corrections to the helium nonrelativistic ground 
energy and the helium ground energy E. E = E + 6EW + 5E^\AE — E — E exp . 
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N 


6E^ 


SE^ 


E 


AE 


20 


-0 00009584 


00002238 


-2 90336 


00003 


30 


-0.00009610 


0.00002239 


-2.90337 


0.00002 


40 


-0.00009586 


0.00002238 


-2.90337 


0.00002 


50 


-0.00009586 


0.00002238 


-2.90338 


0.00001 



It can be seen from Table 1 that the variational ground energy E^ is lower than the 
experimental value and the error will become larger increasing with the number of basis 
functionsiV, which indicates that the Hamiltonian in Eq. (j3J) is not accurate enough. How- 
ever, the case of Eq is opposite to that of Eoo and it approaches to the experimental value 
when N is increased. This shows that the ground energy has been raised after the correction 
of the nuclear motion is considered. One can see from the data in Table 1 and Table 2 that 
it is necessary to consider the relativistic and radiative corrections to the nonrelativistic 
energy. 

Finally, considering all these contributions, we get the high precision value of the helium 
ground energy —2.90338 a.u. and the relative error is 0.00034% (the experimental value is 



-2.90338629 a.u. 



2i|). 



In these calculations, we use a Mathematica software, which is a tool of symbolic calcu- 
lation. All results are exact except for that of the last step that will bring an error. But 
this error can be controlled by setting the precision, and we have ensured that it is much 
smaller than AE. So, the error is mainly derived from the inaccurateness of model, i.e. the 
ignorer of the higher orders relativistic and QED effects. We take the error as half of order 
a 3 relativistic correction and the estimated theoretical uncertainty of calculation is obtained 
to be ±0.00001 a.u.. 

In conclusion, we pointed out that the result of variational calculation of the ground 
state energy for a helium atom would lower than the experimental value when the nuclear 
motion is ignored. First, a new variational wave function with a flexible scaling parameter 
fchas been constructed. Then we take account of the nuclear motion and the relativistic and 
radiative corrections to the nonrelativistic ground energy. Using the variational method and 
perturbation method, the high precision value of the helium ground energy is worked out. 
It is obvious that accuracy will be improved increasing with the number of basis functions 
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N. Thus the ground energy problem of a helium atom is now solved satisfactorily. In our 
later work, higher orders relativistic and QED effects will be considered to agree with the 
observational value better. 

The authors are grateful to Yu-Hong Chen, Yong-Qiang Wang and Zhen-Hua Zhao for 
their useful discussions. 
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